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Abstract. We compare 2 — > 2 covariant transport theory and causal Israel- 
Stewart hydrodynamics in 2+1D longitudinally boost invariant geometry with 
RHIC-like initial conditions and a conformal e = 3p equation of state. The 
pressure evolution in the center of the collision zone and the final differential 
elliptic flow V2 (pt) from the two theories agree remarkably well for a small 
shear viscosity to entropy density ratio rj/s fa l/(47r), and also for a large cross 
section a sa 50 mb. A key to this agreement is keeping all terms in the Israel- 
Stewart equations of motion. Our results indicate promising prospects for the 
applicability of Israel-Stewart dissipative hydrodynamics at RHIC, provided the 
shear viscosity of hot and dense quark-gluon matter is indeed very small for the 
relevant temperatures T ~ 200 - 500 MeV. 



PACS numbers: 24.10.Lx, 24.10.Nz, 25.75.Ld 

1. Introduction 

Recent interest in heavy-ion physics has focused on constraining the transport 
properties of hot and dense nuclear matter using experimental data from RHIC. A 
particular open question is the effect of the conjectured "minimal" shear viscosity 
?7 = s/(4ir) [1] on the dynamics and observables (s is the entropy density). Studies 
of dissipation require a suitable theory framework, in principle a non-equilibrium one. 
However, close to local equilibrium, one can also apply dissipative extensions of ideal 
(Eulcr) hydrodynamics 0G2E]. 

The straightforward Navier-Stokes extension of relativistic ideal hydrodynamics 
with corrections linear in gradients leads to acausal behaviour and instabilities. An 
improved formulation proposed by Mueller and later extended by Israel and Stewart [5] 
(IS) includes second derivatives, which alleviates the causality problem. However, that 
theory originates from an arbitrary truncation of the entropy current at quadratic 
order in dissipative corrections (shear and bulk stress, and heat flow), which is not 
a controlled approximation. Derivations of the IS equations from kinetic theory 
again rely on an arbitrary truncation of nonequilibrium corrections to the phase 
space density at quadratic order in momentum (Grad's 14-moment approximation). 
In contrast, rigorous (Chapman-Enskog) expansion in small gradients near local 
equilibrium results in Navier-Stokes theory. 

Because of these uncertainties about the region of validity of Israel-Stewart theory, 
detailed cross-checks against a nonequilibrium approach are paramount. Here we 
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report on a comparison against covariant transport theory, for conditions expected 
in Au + Au at ^/snn ~ 200 GeV at RHIC, and investigate the effect of small shear 
viscosities on the dynamics and differential elliptic flow V2(pt)- 

2. Covariant transport and Israel-Stewart hydrodynamics 

We solve the equations of motion of causal Israel-Stewart dissipative hydrodynamics 
d^ v = , d^N" = (1) 

Dtt^ = - — (vr^ - 2r)V^u v h - {u^n™ + u v ^)Du a 

- ^(d a u a + Dln^) + 27r^^ x (2) 
T» v = (e + V)vF-u v - pg^ v + tt^ , N» = ra" (3) 

in 2+1D longitudinally boost invariant geometry. Here s, p, u^, and n are the local 
energy density, pressure, flow velocity and particle density; D = u^d^; the () brackets 
denote traceless symmetrization and projection orthogonal to the flow 

A^ v ) = -A^ a A u0 {A afj + A Pa ) - -A^ v A a0 A al3 , A^ = g^ - u^u v , 
2 3 

and V p = A^ v d v . In the equation for the shear stress ir^" we include the term with the 
vorticity tensor ui^ = {\/2)A^ a A u P{dj3U a — d a up) that follows from kinetic theory. 

The numerical IS solutions are obtained with a modified version of the code used 
in [7]. Details of the algorithm will be published elsewhere. Covariant transport 
solutions are calculated via the MPC algorithm [Si [9] . 

To aid comparison with covariant 2^2 transport theory, we take an ideal gas 
equation of state of massless particles, £ = 3p, and keep particle number conserved. For 
our conformal system bulk viscosity vanishes, while heat flow is ignored for simplicity. 
The shear viscosity and shear stress relaxation time of the fluid are matched to the 
values from kinetic theory [3 [6]: n w 6T/(5cr) and tv = 2/3 2 ?7 ~ 9\mfp/§ (because 
/?2 = 3/4p), where an isotropic two-body cross section is utilized for simplicity. We 
consider two scenarios: i) a constant cross section, and ii) a cross section a ~ t 2 / 3 
growing with proper time r = \/t 2 ~ z 2 . During the initial one-dimensional expansion 
stage rj/s ~ t 2 / 3 grows in the former case, while in the latter case r\j s ~ const. 

Initial conditions expected at RHIC for Au + Au at ^snn ~ 200 GeV and impact 
parameter b = 8 fm arc modelled through an initial density profile proportional 
to the (local) number of binary collisions for diffuse gold nuclei, normalized to 
dN/di](b = 0) = 1000 to account for the observed dN c h/dr] w 700. The system is 
assumed to start from local thermal equilibrium {tt^" = 0) at a thermalization time 
To = 0.6 fm with uniform initial temperature To = 385 MeV. 

3. Main results 

Figure [T] compares the evolution of the transverse and longitudinal pressure from 
covariant transport (lines with symbols) and IS hydrodynamics (lines without 
symbols), averaged over the center of the collision zone rr = \/ x 1 + y 2 < 1 fm for 
rj/s rs l/(47r). At such a low shear viscosity, we find that IS hydrodynamics is a good 
approximation to covariant transport in the densest region of the collision. 
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Figure 1. Proper time evolution of the average transverse (T xx , solid) and 
longitudinal pressure (T zz , dashed) in covariant transport (symbols) and causal 
Israel-Stewart hydrodynamics (no symbols) near the center of the collision zone 
(t"t < 1 fin) for RHIC-like initial conditions (see text) and 77/s fa l/(47r). For 
clarity, T zz is divided by a factor of 10. 



In order to calculate momentum distributions, a freezeout prescription is 
necessary. Here we apply Cooper-Frye (sudden) freezeout 

~^-= I da^x)p^f(x,p) , (4) 
dyd z p T J 

where / is the phasespace density corresponding to the local fluid variables and da^ 
is the local normal to the freezeout hypersurface. We choose a constant density 
hypersurface n = 0.365 fm~ 3 , which is the density of an ideal gas of massless gluons 
in chemical equilibrium at T = 120 MeV (3 colours). 

Figure[2jleft) shows elliptic flow from IS hydrodynamics as a function of transverse 
momentum, for the 77/s « l/(4ir) scenario. The dashed curve corresponds to (|4|) 
with the local equilibrium ansatz f eq {x,p) — n{x) / [BitT (xY) exp[— p fJ, u fJ _/T(x)} (we 
use Boltzmann statistics for consistent comparison with transport). This incorporates 
dissipative corrections to the evolution of the flow field, temperature and particle 
number density, but ignores noncquilibrium distortions of the local momentum 
distributions (i.e., assumes -k^ = 0). The solid curve is the result for 

p^p^n^ix) 



f{x,p) = feq(x,p) 



1 



(5) 



8n{x)T 3 (x) 

that properly takes the local shear stress into account. We find that relative to ideal 
hydrodynamics (dotted curve), dissipation reduces elliptic flow by ~ 30%, even at such 
a low 77/s w 1/(4tt). This reinforces an earlier kinetic theory estimate [TU]. Unlike [I], 
we find that at least one third of the reduction comes from corrections to the ideal 
hydrodynamic variables (the lower the px, the larger the fraction), while the remaining 
up to two-thirds come from shear stress altering the local momentum distributions. 

Finally, we compare in Figure fright) differential elliptic flow U2(pt) between 
covariant transport (squares) and IS hydrodynamics (dashed and solid). The dotted 
line is the ideal hydrodynamics reference. For o gg ^g g ~ 47 mb as in [9] we find 
excellent agreement between transport and IS hydro (open squares vs dashed line). 
We also find good agreement for 77/s sa l/(47r) (filled squares vs solid line), for which 
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Figure 2. Left: Elliptic flow as a function of px for Au + Au at 
y/s jvjv ~ 200 GeV and 6 = 8 fm from ideal hydrodynamics (dotted) and IS 
hydrodynamics with Cooper-Frye freezeout ignoring (dashed) or incorporating 
(solid) dissipative corrections to the local momentum distributions (see text) for 
tj/s rj l/(4-7r). Right: Comparison of V2(pt) from covariant transport (squares) 
and IS hydrodynamics (lines) for rj/s l/(47r) (open squares vs dashed), and 
o~ gg — ygg f» 47 nib (filled squares vs solid). The ideal hydro reference is also shown 
(dotted). 



IS hydro somewhat underpredicts the transport results at higher pr > 1.5 GeV. We 
caution that the hydro results at highpx are also sensitive to the freezeout prescription, 
which deserves further investigation. 

The modest ~ 30% dissipative corrections found agree well with [2] and are 
much smaller than the effect in [3]. We suspect that this is because [3] ignored the 
first term in the second line of ([2]). If the system is near global equilibrium, that is 
appropriate. On the other hand, RHIC-like initial conditions have large gradients, and 
our comparison to kinetic theory shows that in that case those terms are important. 

4. Conclusions 

See abstract. 
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